The dearth of halo dwarf galaxies: is there power on short scales? 
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Af-body simulations of structure formation with scale-invariant primordial perturbations show sig- 
nificantly more virialized objects of dwarf-galaxy mass in a typical galactic halo than are observed 
around the Milky Way. We show that the dearth of observed dwarf galaxies could be explained by a 
dramatic downturn in the power spectrum at small distance scales. This suppression of small-scale 
power might also help mitigate the disagreement between cuspy simulated halos and smooth ob- 
served halos, while remaining consistent with Lyman-alpha-forest constraints on small-scale power. 
Such a spectrum could arise in inflationary models with broken scale invariance. 
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The broad-brush picture painted by the inflation- 
inspired hierarchical-clustering paradigm accounts for 
the smoothness of the cosmic microwave background 
(CMB), its tiny temperature fluctuations, the flatness of 
the Universe, and the observed distribution of galaxies. 
However, finer inspection yields some possible — and pos- 
sibly troubling — discrepancies between the models and 
observations that still need to be ironed out. One of 
these is the dearth of substructure in galaxy halos. 

Recent high-resolution A^-body simulations have 
confirmed earlier analytic arguments |^ that suggested 
that hierarchical-clustering models should produce far 
more dwarf galaxies around the Milky Way than are ob- 
served. There are only 11 dwarf galaxies with internal 
velocity dispersions greater than 10 km sec""'^ within the 
virial radius of the Milky Way halo. However, numerical 
simulations of structure formation in a hierarchical model 
indicate that a halo of the Milky Way's mass and circu- 
lar speed should contain roughly an order of magnitude 
more dwarf galaxies. These theoretical results arc robust 
to changes in the values of cosmological parameters or in 
the tilt of the primordial spectrum of perturbations. 

One might at first be tempted to dismiss this discrep- 
ancy between theory and observations as a consequence 
of some nasty astrophysics. After all, the simulations 
consider only gravitational interactions and identify only 
virialized dark halos, while the real Universe is filled with 
gas, and dwarf galaxies are identified by their visible mat- 
ter. So, for example, one might guess that the halos are 
there but remain invisible because the gas has been ex- 
pelled by an early generation of supernovae. However, 
even generous estimates of the efficiency of supernova- 
driven winds fall short of explaining the absence of lu- 
minous dwarf galaxies Further, even if the baryonic 
matter could somehow be driven out of the mini-halos 
or kept dark, we would still have trouble explaining how 
a spiral disk could have formed in the strongly fluctuat- 
ing potential of such a clumpy halo. See Refs. [QJIJ^ for 



more detailed reviews of such arguments. 

In the absence of any prosaic astrophysical mechanism, 
it is natural to think of more exotic explanations. One 
strategy is to modify the nature of the dark matter, in or- 
der to prevent low-mass halos from forming. One option 
is that the dark matter is warm, for example a neutrino of 
mass around a keV, which would suppress the formation 
of small-scale structure by free-streaming out of poten- 
tial wells However, such a neutrino would diminish 
power on scales 2 to 12 h"^ Mpc to a degree that may con- 
flict with the power inferred from the Lyman-alpha for- 
est Another possibility is that the dark-matter 
particles interact strongly with each other, but not with 
ordinary matter [^j. However, the properties required 
of this particle (elastic-scattering cross-sections of order 
10~^^ cm^) are almost inconceivable in the predominant 
paradigms of weakly- interacting massive particles Q and 
axions [ p^ . 

Here we consider an alternative strategy: broken scale 
invariance in the primordial power spectrum, arising from 
some feature in the inflaton potential. In slow-roll infla- 
tion, the amplitude of density perturbations on some co- 
moving scale is proportional to V^/'^ /V , where V is the 
inflaton potential, and V' its derivative, when the comov- 
ing scale under consideration exits the horizon. Suppose 
that V' is initially very small and that there is then a 
discontinuity in the second derivative, V" . The slope V' 
will then jump, and the density-perturbation amplitude 
will drop steeply. In this way, the density-perturbation 
amplitude on some suitably small scale can be suppressed 
relative to the power on some larger scale. Such an idea 
has been invoked as a possible explanation of a claimed 
break in the power spectrum at ~ 100 Mpc scales [p 12 
(and in an attempt produce non-Gaussian features ||13|). 
Here, we use this idea to account for the lack of halo sub- 
structure. Below, we show how a suppression of small- 
scale power can affect the dwarf-galaxy abundance and 
illustrate a particular inflationary model. We also argue 
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briefly that the absence of small-scale power may help 
explain the discrepancy between simulated cuspy halos 
10 and observed smooth halos. 

For simplicity and clarity of presentation, we restrict 
our discussion here to an Einstein-de Sitter (EdS) model, 
but everything can readily be generalized to a low-density 
Universe. Although we focus on an EdS cosmogony, we 
use a cold-dark-matter (CDM) power spectrum with 
r = 0.25 (e.g., as in a rCDM model |16|) to ensure that 
observational constraints from the shape and amplitude 
of the galaxy power spectrum, the cluster number den- 
sity, and the COBE normalization are satisfied 
The solid curve in Fig. |^ shows this power spectrum, 
along with two other power spectra motivated later. The 
corresponding rms mass fluctuations are also shown. 

We need to understand how changes to the power spec- 
trum of density perturbations will affect the abundance 
of dwarf galaxies. We are interested in the abundance 
of mini-halos of some given mass that exist within some 
larger halo of mass 2 x 10^^ Af©, comparable to that of the 
Milky Way and the halo mass used by Moore et al. [Q. 
The simulations of Moore et al. show that the number of 
mini-halos in the Galaxy today is very close to the num- 
ber that existed in the proto-galaxy; i.e., no more than a 
small fraction are fully disrupted during their subsequent 
orbital motion in the Milky Way halo. Thus, our task is 
simplified: we can calculate the abundance of mini-halos 
in the proto-galaxy that later became the Milky Way 
halo. To do so, we use the conditional mass function 
PI, given by 
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This equation gives the fraction of the mass in a galactic 
halo of total mass Mbig that was in mini-halos of mass 
greater than Afsmaii at the redshift Z{ at which the pro- 
togalaxy broke off from the expansion. Here cr(Af) is the 
linear-theory rms fractional mass fluctuation in spheres 
of radii that on average enclose a mass A/, and = 1.7 
is the critical threshold for collapse. 

It is the circular velocity of the dwarf galaxies, not the 
mass, which is observed. We relate the circular speed Vc 
to the mass by assuming the mini-halos underwent col- 
lapse at the protogalaxy-formation epoch. Then the cir- 
cular speed is obtained from = 10 M GH{z[), where 
G is Newton's constant and H{z) is the expansion rate at 
the collapse redshift z. To compare with the numerical 
results of Moore et al. we consider a Galactic halo 
of mass Af = 2 X IO^^Mq with a circular speed of 220 
km sec""'^. Using the conditional mass function, we obtain 
for the CDM power spectrum the cumulative number of 
halos shown by the solid curve in Fig. || (cf. Fig. 1 in 
Ref. [^). Although our calculation of the number of ha- 
los has several shortcomings, the good agreement with 
the numerical results shown in Fig. 1 of Ref. B clearly 
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FIG. 1. The upper panel shows the power spectrum for 
r = 0.25 CDM (solid curve), for a model in which the 
power spectrum is arbitrarily cut off at A; = 4.5 h Mpc (dot- 
ted curve), and the broken-scale- invariance inflation model 
(dashed curve). The lower panel shows the rms mass fluctu- 
ation as a function of the enclosed mean mass M for these 
three models. 



indicates that we are including the essential physics. In 
particular, we reproduce the order-of-magnitude excess 
of halos in the theoretical prediction, as compared to ob- 
servations, for Wc/'^giobai in the range 0.05 to 0.1. 

We now consider how the dwarf-galaxy abundance 
changes if the power spectrum is modified on short 
scales. Clearly, if power is reduced on dwarf-galaxy 
scales (^ 10^° Af©) relative to that at the Galactic scale 
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FIG. 2. The cumulative number of mini-halos for the power 
spectra shown in Fig. |l| as a function of the circular speed Vc of 
the halo divided by the circular speed Ugiobai of the Galactic 
halo. The points show the Milky Way satellites. Compare 
with Fig. 1 in Ref. g. 

10^^ Mq), the number of dwarf galaxies will be sup- 
pressed. To investigate what is required, we first sharply 
cut off the power spectrum at fc = 4.5/i Mpc~^, shown by 
the dotted curve in Fig. the lower panel shows how the 
rms mass fluctuations are modified. The corresponding 
prediction for halo substructure is shown as the dotted 
curve in Fig. ^ and is in striking agreement with the 
observations. 

However, such a cutoff is purely phenomenological, and 
much stronger motivation is needed. In particular, the 
question arises as to whether inflation can produce such 
a power spectrum. The heuristic arguments given above 
suggest that the power spectrum can undergo a sharp 
drop on small scales if there is a discontinuity (or near- 
discontinuity) in the slope of the potential, and this is 
known as the broken-scale-invariance (BSI) model. Re- 
markably, there is an exact solution for the power spec- 
trum in this situation, due to Starobinsky ||20[| , even 
though the slow-roll approximation does not apply. Fur- 
ther, this spectrum has a universal form, depending only 
on the change in the derivative across the discontinuity, 
in the sense that the form of the spectrum is preserved if 
the discontinuity is smoothed out, provided the smooth- 
ing is over a sufficiently short range of scalar-field values. 
The underlying reason for this is that the field briefly be- 
comes kinetic-energy dominated as it goes over the dis- 
continuity, and the precise form of the potential becomes 
irrelevant. 

The dashed curve in Fig. |l] shows this power spectrum 
for a suitably chosen location and amplitude of the dis- 



continuity (the amplitude is p = 10, in the notation of 
Ref. The universal form features a modest rise, 

followed by a series of oscillations asymptoting to the 
original spectral shape at a much lower amplitude; in 
this case the power is reduced by a factor — 100 on 
short scales. Such a model represents the fastest possible 
cutoff in power that can be obtained from a single-field 
inflation modelj^ 

This power spectrum gives rise to the dashed curve 
shown in Fig. |[ This model is almost as successful as 
the cutoff power spectrum; it produces roughly 20 mini- 
halos with Wc/vgiobai ^ 0.04, in much better agreement 
than the original power spectrum with the 0(10) that 
are observed. We conclude that the BSI model provides a 
promising possibility for reconciling the predictions with 
observations. 

Our modification to the power spectrum sets in only 
on very short scales, and so does not affect successes of 
the standard paradigm on much larger scales, such as 
the cluster number density or galaxy correlation func- 
tion. However, we need to check that the loss of short- 
scale power is not inconsistent with object abundances at 
high redshift. The most powerful and direct constraint 
comes from the abundance of damped Lyman-alpha sys- 
tems |2^. Observations exist at redshifts 3 and 4, giv- 
ing comparable constraints. Redoing the calculation of 
Ref. ||l^ to include updated observations ||2^ gives a 95% 
confidence lower limit of 

ct(1O^"/i~^M0,z = O) > 2.75 + /1, (2) 

in a critical-density Universe, which is (just) satisfied by 
the BSI model (though not by the cutoff power spec- 
trum). 

This narrow escape makes it look as if the model is 
quite marginal. However, we have only analyzed the 
critical-density case, and the damped-Lyman-alpha con- 
straint is much weaker in a low-density model. The ob- 
servational constraint on a{M) at redshift 4 is almost the 
same (no more than 10% lower for f2o values of interest 
p3| ) . However at redshift 4 suppression of the growth of 
perturbations has yet to set in, and the normalization to 
COBE (and/or cluster abundance) means a is higher by 
a factor l/fto- With this additional factor all our models 
would easily pass the damped Lyman alpha system test 
were flo ~ 0.3. 

There are two other key short-scale predictions worth 
considering. One is the power spectrum of perturbations 
obtained from the Lyman-alpha forest However, 
the power spectrum derived in Ref. H does not extend 



*We suspect that this is true in multi-field models too, as 
additional fields will supply extra friction, slowing down the 
fields' evolution and hence stretching features in fc-space. 
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much beyond k = 5h Mpc""'^, which is where our model 
begins to differ from the usual one. It therefore does 
not appear in conflict with the power in our theoreti- 
cal spectra; indeed, if anything, their analyses see less 
power than expected on short scales. The second test is 
whether the model has early enough structure formation 
to reionize the Universe by a redshift 5, as required by 
the Gunn-Peterson test. Estimates using the technique 
of Ref. j25| would suggest this is marginal in the critical- 
density case. However, once again in the low-density case 
we benefit from the much higher early-time normalization 
of the power spectrum, which increases the predicted red- 
shift by a factor ~ l/Oo- We thus conclude that this BSI 
model is likely to be consistent in the observationally- 
favored low-density flat model, though a more careful 
investigation is merited. 

The suppression of small-scale power may also help to 
explain the discrepancy between the cuspy inner halos 
observed in simulations and the smooth halos inferred 
from observations. It has been argued that the steepness 
of the halo density proflle at small radii depends on how 
the characteristic density of the mini-halos that merge to 
form galactic halos scales with their mass. If small halos 
are denser — as they would be if they collapsed earlier — 
then the cusp is steeper p^ . If power is suppressed on 
sub-galactic scales, then all mini-halos with sub-galactic 
masses will undergo collapse at about the same time and 
thus have similar densities |^^. If so, then halos should 
be less cuspy than prior simulations have shown. How- 
ever, these heuristic arguments have been tested numeri- 
cally for the case of galaxy clusters , where no change 
to the core structure was found when short-scale power 
was cut off. More work will be needed assess the effects 
on galactic halos. 

In conclusion, the arguments |ll|-^,H suggesting that 
the standard structure-formation paradigm predicts too 
much substructure in galactic halos have now been re- 
fined to an extent that requires that the problem be taken 
very seriously and that radical answers to the discrepancy 
be considered. We have proposed here that the discrep- 
ancy indicates a dramatic lack of short-scale power in 
the primordial power spectrum. While such a sharp fea- 
ture in the spectrum would not be expected a priori, it 
can be obtained from infiation, with the BSI model be- 
ing the best available example. This model has already 
been invoked to introduce features on much larger scales 
of around 100/i~^ Mpc [pl] , p^ , but the observational 
motivation for a short-scale break appears stronger. If 
our suggestion is correct, the lack of short-scale power 
will soon become apparent in new observations of high- 
redshift phenomena. 
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